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I. INTRODUCTION 



The problcm of the classical limit of quantum mechanics has a long history. In the beginning, on the basis of the 
analogy with special relativity where the limit c — > oo leads to the classical behavior, it was thought that the classical 
limit was just the limit H — > 0. But it was soon realized that this was only one element of the problem, namely, 
macroscopicity, and that other elements must be taken into account: e.g. quantum mechanics has a probabilistic 
non-Boolean structure while classical mechanics has a non-probabilistic and Boolean one. Thus, necessarily two new 
i | elements must come into play: 

C ■ 

ç$ • Decoherence, that transforms the non-Boolean structure into a Boolean one, and 

3 ! 

^jH, • Localization (actualization or the choice of a trajectory) that, with macroscopicity -which circumvents the un- 

r* ■ certainty principle- turns the probabilistic structure into a non-probabilistic one. 

• i— i , 

In general, decoherence in quantum systems is defined as a process that leads to the diagonalization of a density 
matrix (more precisely, to the decay of the cross-terms in the expectation value of an observable in some basis). In a 
first period, decoherence was explained as the result of the destructive interference of the off-diagonal elements of the 
density matrix (see [1], [2]); however, this line of research was abandoned due to technical difíiculties derived from 
the formalism used to describe the process. As a consequence, decoherence begun to be conceived as produced by 
the interaction between a system and its environment. This approach gave rise to the einselection program, based 
on the works of Zeh ( [3], [4], [5]) and later developed by Zurek and coworkers ( [6], [7], [8], [9], [10], [11], [12]). 
Although many relevant rcsults have becn obtaincd by means of einselection, this approach still involves certain 
unsolved problems, as those related with the explanation of the emergence of classicality in closed quantum systems, 
the criterion for introducing the 'cut' between the system and its environment, and the defmition of the preferred 
('pointer') basis where the system behaves classically (see [13], [14]). As the result of these and other difíiculties, a 
number of alternative accounts of decoherence have been proposed (see [15], [16], [17], [18], [19], [20]). 

In a series of papers ( [21], [22], [23], [24], [25], [26], [27], [28], [29], [30], [31], [32], [33], [34], [13]) we have returned 
to the initial idea of the destructive interference of the off-diagonal terms of the density matrix, but now on the basis 
of a different formalism: the formalism introduced by van Hove ( [36], [37], [38], [39], [40]). We have called this new 
approach ' self-induced decoherence' [13] because, from this viewpoint, decoherence is not produced by the interaction 
between a system and its environment, but results from the own dynamics of the whole quantum system governed 
by a Hamiltonian with continuous spectrum. In this approach, the difíiculties derived from the einselection program 
are absent: self-induced decoherence can be used in closed systems as the universe [30], the defmition of a convenient 
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subalgebra plays the role of the coarse-graining induced by the environment, avoiding the 'cut' problem [13], and the 
pointer basis is perfectly defined (see [13], [34] and Section III. C below). 

Sclf-induced decoherence is capable of addressing relevant problems from a general perspective, e.g. the problem of 
supplying a good definition of the classical limit in all cases 1 . Let us explain the essence of the idea supporting this new 
approach. When we deal with continuous spectra, the destructivc interference is embodied in the Ricmann-Lebesgue 
theorem which states that, if f(y) G Li, then 

lim / f(v)e- VJ idv = 

where e~ li # is the v— oscillating factor that produces the destructive interference. In the case of decoherence, v = 
lü — u' , where uj,lo' are the continuous indices of the density operator p(u),u)') in the energy eigenbasis; then, v = 
corresponds to the diagonal. However, to require that f(v) G Li is to ask too much regularity to function f(v), 
because in this case not only the off-diagonal [y ^ 0) terms, but also the diagonal [y = 0) ones will vanish when 
t — > oo. Thcrcforc, we use less regular functions, precisely f(v) = A8{v) + where only fi(v) G Li. In this case, 

/a 
f{v)e~ %! %dv = A 
-CL 

and the diagonal terms A8{v) remain while the off-diagonal terms f\{v) vanish. We will apply this main idea to 
the non-integrable case, and present the computations in all detail in Section III. B, by using our previous results on 
quantum systems with continuous spectrum containcd in papers [22], [23], [24], [27], [28] and [29]. With this strategy 
we have already obtained, in paper [26], the classical limit for integrable systems. We have also presented this result 
in more rigorous mathematical basis in [33] and cxplaincd the physical foundations of the method in papers [13] and 
[34]. The mathematical basis of the theory is explained in papers [41] and [32]. But, of course, the big challenge to 
prové the consisteney and generality of the method is to find its version for non-integrable systems, obtaining unstable, 
potentially chaotic classical trajectòries as a final result, which could explain models as those of ref. [42]. 

In the case of integrable systems, the classical limit was obtained by a combination of the van Hove formalism and 
the Weyl-Wigner-Moyal isomorphism in a globally defined pointer basis. But in the non-integrable case, such a global 
basis does not exist. Nevertheless, the just quoted isomorphism is what allows us to relax the global condition and 
to generalize the formalism: quantum mechanics is formulated in a phase space that is covered with charts where 
local pointer bases can be defined. The set of all these local pointer bases will yield decomposition (38), which is the 
essential tool of this paper. 

The formalism of the theory is presented in a self-comprehensive way, with a mathematics as simple as possible and 
in the simplest possible case; this seems enough for the physical purposes of this paper. In Section II, a brief review 
of the Weyl-Wigner-Moyal mapping is developed, and in Section III, the theory of decoherence in non-integrable 
systems is explained. In Section IV, the classical limit of quantum non-integrable system is obtained. In Section V, 
the localization phenomena is bricfly discussed, and in Section VI, our previous results are generalized to the case 
of partially non-integrable systems. In the conclusion (Section VII), we list the possible future applications of the 
theory and explain why it could be considered as a minimal formalism for quantum chaos. Finally, in Appendix A 
we explain the integration of two systems of diffcrential equations relevant to our theory, and in Appendix B we give 
an example of non-integrable system. 

II. WEYL-WIGNER-MOYAL MAPPING 

Let M. = M 2 (n+i) = K 2 ' Ar+1 ^ be the phase space of our classical system. The functions over this phase space will 
be called f(<j>), where <j> symbolizes the coordinates over Ai 

<j> a = (q 1 ,...,q N+1 ,p 1 q ,...,pï +1 ) o=l,2,...2(JV+l) (1) 

As it is known (see [43], [44]), we can map A, the àlgebra of regular operators O of our quantum system, on A q , the 
àlgebra of integrable functions over M , via the Wigner symbol 



1 Precisely: in all cases where the system do have a classical limit; e.g., systems with no quasi-continuous limit, yielding to a 
non-continuous energy spectrum, are excluded [35]. 
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symb : À — ► A q symbO = 0{<j>) (2) 
Precisely: let us consider that M has a symplectic form 

Then, 

symbf 



= /(0) = J d 2(W+1) V> exp (j^°üW 6 ) (f (V)/) (4) 



where / e À, f{4>) E A q , and 



i 



T(V) = cxp ^°w o6 6 J (5) 

On we can definc the star product (i.c. the classical operator related with the multiplication on A and, therefore, 
dcfining the corresponding operation on A q ) as 

symb(fg) = symbf * symbg = (/ * g)(4>) (6) 
It can be proved ( [43], eq.(2.59)) that 

(/ * 9){<t>) = f{4>) exp (-yX^ ab ^i) .9(0) - 5(0) exp (^*d w o6 "2 6 ) /(0) (7) 

We also define the Moyal bracket as the symbol corresponding to the commutator in A 

{f,g} mb = ±(f*g-g*f) = S ym6 (^[/, <?]) - 7^/(0) sin ("f V^) .9(0) (8) 

In the limit fi — > 0, the star product becomes the ordinary product, and the Moyal bracket becomes the Poisson 
bracket 2 

(/*<7)(0) = /(0M0) + O(fi) (9) 

{/,sW = {/,s} P fc + o(fi 2 ) (io) 

Then, we can either say that when fi — > the quantum structure tends to the classical one, or that when fi ^ the 
classical structure is quantized or deformed into the quantum one. 
Let us observe that if / commutes with g, eqs.(7) and (9) change to 

(/ * .9X0) = f{4>) cos (-yKw" 4 ^) .9(0) (11) 
(/*ff)(0)-/(0).9(0)+O(fi 2 ) (12) 

as it can be proved using eq.(7). 

Finally, if we want that the mapping symb be one-to-one, we must define a unique inverse of symb, namely, the 
usual quantization rule q — > q , p — > p cndowed with a unique ordering prescription, e.g. the symmetrical or FFej/Z 
ordering prescription that maps 



2 From cq. (7) it is clear that these 0(fi) are continuous functions in the limit h = 0. This fact will be important in Section V. 
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symb 1 {qp) = - (qp + pq) 



(13) 



Then, we have 

symb -1 : A q -^>A, symb : A-^> A q (14) 

The one-to-one mapping so defined is the Weyl-Wigner-Moyal symbol. With symb -1 ' we can 'deform' the classical 
system and obtain a quantum mcchanical systcm. With symb we go from usual quantum mechanics to a quantum 
mechanics 'alia clàssica', formulated over a phase space Ai, that will become the usual classical picture in the limit 
h — > (as we will explain below in detail). The relation between the two structures, given by eq.(14) (and eq.(66) 
below), is an isomorphism that we will call Weyl-Wigner-Moyal isomorphism, the only one we will use in this paper. 

Since A is a space of operators on a Hilbert space H, so it is its dual A'; then, as it is known, the symbol for any 
p E A' is defined as 3 

p(<f>) = symbp = (2tt hy {N+1) symb (íor op orator S )P (15) 
where the symb for operators is defined by the eqs.(4) and (5). From this definition, we have (see [43], eq.(2.13)) 

(p|0) = (symbp\symbd) = j d(f> 2{N+1) p(4>)0((f>) (16) 

and in A and A' all the equations are the usual ones (i.e. those of papers [43] and [44]). Let us remark that the last 
cquation is the cornerstone of our theory of the classical limit. In fact, as we will see, it will remain the same when we 
go from regular to singular objeets. Once this statement is understood, the translation from the quantum language to 
the classical one will be easy. 



III. DECOHERENCE IN NON-INTEGRABLE SYSTEMS 



A. Local CSCO 



a.- We will begin with demonstrating an important theorem: when our quantum system is endowed with a CSCO 
of TV + 1 observables containing H that defines a basis in terms of which the state of the system can be expressed, the 
underlying classical system is integrable. In fact, let a classical system be defined in a phase space Ai =M 2 ' A ' +1 - ) that 
can be deformed 'alia WeyP. If our quantum system is endowed with a N + 1— CSCO {H,Oi, ...,Ò n}, the Moyal 
brackets of these quantities are 

{O 7 (0),O. 7 (0)} m6 = symb (^[Òj^j^j = (17) 

where /, J, ... = 0, 1, N and H = Oo- Then, when h — > 0, from eq.(10) we know that 

{Oj(0),Oj(0)} p6 = O (18) 

Thus, as H(<f>) = Oq (</)), the set {Oj(4>)} is a complete set of N + 1 constants of the motion in involution, globally 
defined over all Ai and, therefore, the system is integrable, q. e. d. 

As a consequence, non- integrable classical systems (precisely classical and also macroscopic ones such that fi«0), 
in their quantum version, cannot have a CSCO of N + 1 observables globally defined containing H. But, according 
to the self-induced approach, the pointer basis is precisely the eigenbasis a global N + 1— CSCO (in such a way 
that the vectors of the pointer basis turn out to be stationary states, see [26]). Therefore, pointer bases cannot be 



3 In the case of states, we must add a new factor (27r/ï) _ ' iv+1 · ) to definition (4) in order to preserve the usual normalization of 
p(4>). However, p{4>) is not non-negatively defined. With decoherence and Ji^Owe will obtain a non-negatively defined p(4>), 
and Aq — > A, the classical boolean àlgebra of Li operators over M. 
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globally defined in non-integrable systems. These systems can be adequately quantized, but it is impossible (at least 
globally) to define a complete stationary eigenbasis of N + 1 — CSCO and, a fortiori, a pointer N + 1 — CSCO or a 
pointer basis where the system would decohere according to the self-induced approach. 4 This is the main problem 
with non-integrable quantum systems. 5 

b.- We will now prové that N + 1 constants of the motion in involution always exist locally. 6 Let us consider a 
non-integrable quantum system (i.e. with no global N+ 1 — CSCO), but let us suppose that, as usual, H{<j>) — symbH 
is globally defined over Ai (this means that any non-global CSCO has at least one global observable: H). 7 Now we 
can try to find N constants of the motion {Oi(<f>)} (7=1,2,..., N) satisfying 



This is a system of N partial diffcrcntial cquations which, with adequate boundary conditions, has a unique solution 
in a maximal domain of integration T> ( j >i around any point <f>i £ Ai (provided that the functions involved satisfy 
reasonable -e.g. Lipschitzian- mathcmatical conditions that we assume 8 ). 

But we would like to obtain a set of constants of the motion in involution. Then, let us suppose that N different 
initial conditions for eq.(19) are given in a 2N + 1 dimensional hypersurface containing fc, that we will call T>^.. 
Integrating (19) we will obtain N constants of the motion Oi(cj>). Moreover, we can easily show that, if these solutions 
are in involution in £>í , they will remain in involution in the domain V t j >i = T>¥ +1 of 2{N + 1) dimensions. In fact, 



Then, since Oi{4>) and Oj{4>) are constants of the motion in P^, the {Oi (</>), Oj((j))} p b will also be so. As a 
consequence, if we could define N constants of the motion such that 



of dimension 2(N +1), as promised. 

Now the problem is reduced to prové the existence of the N Oi((j)), Oj{<p) satisfying eq.(21) in V^.. Again, the 
existence of such a set can be easily proved by using the same strategy as above, but now recursively. We can begin 
with an arbitrary function 0\{4>) defined in a domain D9. of N + 1 dimensions. Then, we consider another 02(4>) 
(defined in a N + 2 dimension domain T>^, containing T>^.) as the Hamiltonian of eq.(19) and obtain by integration 
a function Oi (</>), defined in the domain V^. of N + 2 dimensions, such that in this domain {02(<j>),Oi(<j>)} p b = 0. 
Finally, we iterate the procedure up to find the set of functions in involution in the T>^. of dimensions 2N + 1, which 
can be taken as initial conditions of eq.(19). In this way, the proof is completed. 

c- Now, in order to go from classical to quantum, we can also extend these local Oi(<j>), defined in — 
of dimensions 2(N + 1), to all Ai by defining Oi(4>) = for G AÍW^^ In this case, there will be a jump in the 
frontier of V^, and the definition will be only continuous a.e. (almost everywhere). Or, on physical grounds, we can 
take the precaution of joining these zero functions with functions Oj(</>) in a zone around 2?^, that we will call J 7 ^, 
in an smooth way (e.g. by using C r functions with an adequate r). 

Therefore, we have proved the existence of local complete systems of constants of the motion in involution 
{Oi{<p)} = {H{4>), Oi (0), On{4>)} that we can extend to all Ai, at least a.e., by adding null functions in Ai\ 
2>^ +1 as explained above. Since they belong to 2?^, we will call them {H{(j)),0 ( j )i i{4>),...,0^ i N{<Í>)}· Each system 
{H ((/)), 0,^1(0), O c f >i ^{(j))} can be considered as a local (approximate) N + 1— CSCO in V ( p i = 2?^ +1 in the sense 
that, even if it is not an exact CSCO, we can compute their Weyl transformations obtaining 



4 Observe that, if the CSCO has < N + 1 operators, we have not good quantum numbers enough to label the eigenvectors. 
5 In the 'old quantization' approach, the problems were certainly more severe. 

6 This fact can be considered as almost evident, but since it is not demonstrated in usual textbooks, we will give a complete 
demonstration below. 
7 We can also say that the quantum system is dissymmetrized [45] . 

8 See Appendix A for details. Moreover, a certain determinant A, defined in [46], must be A 7^ in this domain. 




(19) 




(20) 
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{HfcjOfci, O^n} 

and their Wigncr transformations are a complete set of constants of the motion in involution in P^. In fact, from 
eq.(10) we see that 

{O4M, 0^.j{<P)} mb = 0(h 2 ), or [0^,0^ j] = 0(h 2 ) (22) 

namely, they only commute approximately. 

Let us now consider in more detail the joining zones where we have used C r -functions that do not satisfy the 
required differential equations (19) to (21), in such a way that the terms ^ d a oj ab d b of eq.(8) produce unwanted 
contributions of order h/PQ, where P and Q are of the order of magnitude of the jumps in the momentum and 
configuration variables in the joining zone. Since PQ — e 2 is an action measuring the joining zone (where e is the 
characteristic mean width of the joining zone, precisely e 2 ( 7V+1 ) = V e , the volume of the joining zones J^J, the 

unwanted terms are of the order of h/e 2 , that is, they are another contribution 0(|r), or simply 0(ft 2 ), to add to (22) 9 . 
Anyhow, thcse terms will vanish when we make the limit h — > in Section IV. 

Although this approximation seems sufRcient for physical purposes, we can even improve it. In fact, we can repeat 
all the process based on eqs.(19) and (21) substituting them with 

{H(cf>), 0!(cf>)} mb = 0, {0/(0), Oj(cf>)} mb = (23) 

Taking into account eqs.(7) and (10), these are differential equations of infmite order (see Appendix A.c). But we can 
cut these expansions at any finite order and obtain a system of usual finite differential equations that can be solved; 
then, we can repeat the procedure described above with any desired precision. In this way, we can eliminate the 0(h 2 ) 
coming from cq.(10) but not those coming from the joining zones. Thus, we have again obtained approximate (up to 
0(h 2 )) local N + 1— CSCO, 10 and we can define local eigenstates and write equations like (35) for each X*^. 11 

d. - Lct us observe that natural global coordinates <p = (q>Pp) 01 phase space M can be (locally) substituted, by 
using (local) canonical transformation, with (local) coordinates (0..i, O^i), with i = 0, 1, ...N and H — O^o, where 
the 0,,i((j>) are the coordinates canonically conjugated to the O^j^) in 2?^. The (0,.i, 0<pj) is clearly a chart of M. 
in the domain 2?^. 12 Since the system is endowed with adequate smooth properties (let us say C r ), another similarly 
constructed chart (O^ji, 0<f>ji) in the domain 2?^ is smoothly connected with the previous one at any <j) G 2?,^ H V^ >j 
(see demonstration in Section V). Then, the set of all these charts is a C r — atlas in M. This will be the atlas we will 
primarily concerned with. 

e. - We can also dcfinc a (ad hoc) positive partition of the identity (see [51] sec. 3.4) in the following sense. Let us 
define 



9 Any counterfactual 0(h) is really a factual 0(h/S), where the action S — > 00. Precisely, this means that S >> So, where S'ò 
is a characteristic action. In order to give an example for this So, we can consider the fine structure constant: 

h h r _ 1 

m e ca So 137 

where m e is the electron mass, c is the velocity of light, and ao is the Bohr radius. For 5 > 5o we can neglect a and, therefore, 
we loose pure quantum efïects, like spin, and pass from the realm of quantum mechanics to the classical limit. 
10 Counterfactually, O(^) goes to zero when H — > 0. Factually, it goes to zero if the action of the system S is infinitely large. 

Then, since e 2 has the only constraint e 2 < S, it can be as large as we wish and we have ^ — > (see also footnote 13). 

n An example of this phenomenon is the Sinai billiard discussed in Appendix B. Other examples are classical scattering 
systems: in fact, they have an 'in' CSCO and an 'out' CSCO, which are different since the constants of the motion are not the 
same in these CSCOs. Another example is the two slits experiment when we mimic the screen with an infmite potential wall: 
before the screen we have a local CSCO {H, P}, on the screen the CSCO is {H, X} since the kinetic term of the Hamiltonian 
can be neglected with respect to the infinite potential wall, and after the screen again {H, P}. More complex examples are the 
so-called pseudointegrable systems ( [47], p.98, [48], [49]). Tori become spheres with 'handles' that cannot be covered with a 
single chart. A further example is Robnik's billiard [50]. Moreover, it is clear that the fractal structure of some examples of 
chaos breaks the tori completely and, therefore, in this case the radius of the integration domains probably vanish (nevertheless, 
this would not be a physical case, see Appendix A). 

12 This is not a genèric chart, but a very peculiar one, since coordinates O^/ are constants of the motion satisfying eqs.(19) 
and (21). 
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l = I(0) = Çl^(0 (24) 

i 

where B^ >i (<ft) are 'bump' functions such that 

f = i if <f> e 

B* 4 (^H e [0, 1] if € F 0< (25) 
l = if £ U F 0í 

where D ( p i is a domain and F ( f )i is the frontier zone around D t j >i (the are similar to the T^ >i but they are related 
to the -D^J defined in such a way that D ( p i U F < f >i C and the intersection zones of the D's vanish: D ( f >i n -D^ = 0. 
Let us stress that the B^ i {<j>) in the frontier zones satisfy eq.(24). Now, for any A((f>) we can define a 

and for any A{<j>) we have 

A(<f>) = AMY^BtM = Ç^(0) 

i i 

With the mapping symb^ 1 wc find 

À = aymb- 1 A(<f>) = J2 symb^A^ (0) = ^ (26) 

where — symb~ l can be considered as a localization of A in Z)^ 4 . Then, from eq.(26) 

a = J2 à ^ ( 2? ) 

i 

Moreover, since we have a local N + 1— CSCO in each D ( p i U i 7 ^ C X>^, we can decompose 

^E^tf^OI^ (28) 
i 

where the {|i)^} are the corresponding eigenvectors of A^; the j can be considered as a local N + 1— CSCO 
oi' />,, 7\. . 

Now we can prové that the support of symb\j)^ is contained in D^ >i U F < p i i.e. the support of symbA^. In 

fact, from eq.(28) we have 

= A j<t>i \j)^ 

or 

Then, 

S ymòA 0i * symò|j)$f (j^ = A^sym&lj}^ (29) 

But symbA c f >i and all its derivatives vanish for <f> D ( p i U F < f >i . Therefore, if Aj^ 7^ 0, this also must happen for 
symb\j)^ and the support of this function is contained in D ( p i Uf^. If Aj^ = 0, we can repeat the argument 

with the operator A ( p i + aB^ >i and take the limit a — > 0, and we will find the same result. 
From eq.(28) we have 

^ = E^b·)if01í ) (30) 

ij 

all over M. Moreover, from eq.(28) we also have 
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symbÀ^ =Y^ A iïi s y mb \j)^(j\ { <h ) 
i 

and, as we have just proved, 

symb^ij]^^) = if ^ I>* C Dfc U Ffc 
Thcn, since for i =/= fc, D<p z D^, k — (but F ( p i n F ( j >j ^ 0), we have 

ioiííWi 8 = o·i!í ) ii'>ffo·'iií ) ii>5? ) - (liO^o-'iifib^foiif ) = 

/ ^li^OliíWlj^OlSíV^ = 0(e 2(AÍ+1) ) (31) 

where F is the union of all the joining zones F t j )i and e is the characteristic width of the joining zone. Therefore, for 
i^k and e — > 0, 13 we obtain 

OlifV^O (32) 

This means that, in the limit e — > 0, decomposition (30) is an orthogonal decomposition in the Ij}^ 4 ''- 

Remark. Let us now discuss the subspaces that can be defined by the above decomposition. Bcing {\j)^} a 
basis of the Hilbcrt space TL where we are working, any e Ti could be decomposed as 

\<P) = ÇvíjIÍ)^ = Çbi) 

where 

and this ket belongs to subspace Hi oïH. Then, 

W=0H, (33) 

i 

Neverthclcss, from eq.(28) where the off-diagonal terms i ^ j arc absent, we have 



13 Precisely: let us call Vm the volume of phase space: Vm ~ S N+1 . Analogously, 

I M = J^ai^bO^Ol^WW^OI^W 2 ^ - Km ~ 

Let us also define I £ — J T symb\j'}^ (j'\^symb\ j}^ d(f> 2< · N+1 ï ~ V e = e 2 ^ AÍ+1 \ In ordcr to prové eq.(31), it is necessary 

that J £ <C /ai in such a way that / e could be neglected. But I e ~ V e and /ai ~ S' JV+1 ; thcn, e 2 « S . 
Therefore, e must be: 

1. - Such that the ratio \ be negligible to eliminate the unwanted terms ^ *d a u ab ~eí b in the joining zone (see footnote 10). 

2. - As small as e 2 < S to satisfy eq.(31). 

Since h -C S, we can satisfy both conditions with an adequate e, namely, such that: 

« e 2 < 5 
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A+, g n- L ® H t = o t 



and from À and eq.(27) we have 14 

= ^0^=0^ = 0^®^ (34) 

This shows that there are no cross terms Hi <8> Hj in the decomposition of O = H <g> W. 

We can see that the decomposition that really matters for our discussion is (34), the decomposition in subspaces Oi, 

and not (33): the repeated index i in the basis means that this basis corresponds to the decomposition 

donc in the Oi, which is the relevant one for this paper. 



B. Decoherence in the energy 

a. - We will now introduce decoherence according to the self-induced approach. Let us dchnc, in each D^, a local 
N + 1— CSCO where, as in eq.(30), the observables of the N + 1 — CSCO {H, O^} are decomposed as 

H= I u>y^ j \üJ,m) (j , i {u),rn\ l i, i <L·j, O^j = / y^Q TOf ^,. \uj,m) í f )i {LJ,m\ l f )i duj (35) 

where the energy spectrum is < oj < oo and mj^ >i = {mi^, mA^}, mj^ 6 N (the spectra of the O^i are discrete 
for simplicity). 15 Therefore 

Hlu)^)^ = oj\u>, m}^, Ofcilu,™,)^ = O mi ^, (Wjm)^ (36) 

where the \tú,m) l j )i are the eigenvectors of the observables H, and O t j >i (such that symbO^ = 0^(0) ^ only in 
Dcpi HF^J and m is a shorthand for m^j — {m^i, to^jv}. The set {\u> , m) ^ is orthonormal in u) and in m, in 
the usual eigenvalue indices and in i, as proved in eq.(32): 

(uj,m\ ( j >i \u)',m') < j >:j = S(üj - u}')5 mm >5ij (37) 

b. - Now we can dehne our relevant àlgebra of observables. This choice will play the role of coarse-graining in our 
approach. A genèric observable reads, in the orthonormal basis just defined, 



/>oo />oo 

O = ^ / / dL·jdL·j'0(iü,L·j')^ imm/ \L·J,m) l f )i {Li)',m'\ l f )i 
, mm > Jo Jo 



(38) 



where 0(Lú,Lü')<p imm i is a genèric kernel or distribution in lo, ui'. 16 But we must restrict this set of observables since 
it is too large for our purposes; furthermore, it is not easy to work with genèric kernels or distributions. However, 
we cannot make the àlgebra too small either. In fact, let us suppose that, in order to make computation easier, we 
postulate that the 0(u>,uj')^ imm ' be just regular functions. Then, the states read 

poc p oo 

P= y2 / / dw^'p(w,Cj')0, m m'| W : m )0(^': m 'Ui 

, Jo Jo 

mim' 

where the p(w, w')^ mm ' , in the dual space, are also regular functions. Then, 



14 This decomposition is similar to the decomposition of a function in its even and odd parts in a Fourier transformation, 
where the sin is the even basis and the cos is the odd basis. 

15 Hamiltonians with continuous spectra are considered in papers [23] and [24]. We use this kind of spectra since they are 
the usual ones in the macroscopic limit h — > (see [47] eq.(3.1.24) p.67). Strictly, we should call \iü,m)^ i ( ír ·°«»' the vectors 
\iü,m) r f >i , but we will just call them \uj,m)^ i for simplicity. 

As explained at the end of the last subsection, the index i in projector \ui, m)$ i {lo', m!\$ i corresponds to the fact that the 
decomposition is done in the Oi and, therefore, the index is repeated in \tü,m) r f >i and in (o/, m']^ . 
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and, since the product p{ui, oj')fa mm '0(u>, to')famm' is a regular function (i.c. Li in v = ui — lo'), as a result of the 
Riemann-Lebesgue theorem the mean value {0)^u\ would vanish for t — > oo: we would obtain destructive interference 

not only for the off-diagonal terms, but for all of them. On the contrary, if p(u), w')^, imm ' and 0(u>, w')^ mm < were 
genèric kernels, we could not use the Riemann-Lebesgue theorem, and we can presume that there will be no destructive 
interference. This means that 0(u>,u>')fa mm i cannot be so regular nor so non-regular: we must choose something in 
between. In order to avoid these unacceptable results, the simplest choice is the van Hove choice; so, as in paper [26], 
we will takc: 

0(uj, u')fa mm > = 0(u>)fa mm >S(u> - uj') + 0(üj, u')fa mm > (39) 

where the 0(u),u)')fa mm i are ordinary functions of the real variables lo and lo' (these functions must have some 
mathematieal properties in order to devclop the theory; these properties are listed in paper [24]). This choice is 
theoretically explained in papers [36], [37], [38], [39], [40], [41], and [26]. Moreover, we need the 5{lü — üj') term in 
order that the members of the N + 1 — CSCO of eq.(35) be contained in the space of observables. So our operator 
belongs to an àlgebra A (defined by eq.(39) and the properties just required for the O(oj,uj')fa mm '), and reads 



poc poc poc 

O = ^2 / dwO(uj)fa mm '\Li),m)fa(uj,m'\fa + V] / / a^a^'0(u},w')fa mm i\w,rn)fa(uj' ' ,m'\fa 



(40) 



The first term in the r.h.s. will be callcd Os, the singular component, and the second term will be called Or, the 
regular component, 17 and [H , Os] = 0. The observables are the self-adjoint 0^=0 operators. We will say that these 
observables bclong to a space O (which is contained in the operator àlgebra A ); {\ui, m, m')^, \oj,oj' ,m,m') ( f )i } is a 
basis of this space, where 

\uj,m,m') t j >i = \w, Tïij^uj, m'ifc, \uj,J,m,m')^ i = \uj,m) lt)i {u) l ,m'\ <l , i (41) 

Then, the classical analogue of eq.(40) would be 

poc poc poc 

0{4>) = I d^O^famm' \u), m, m' (</>)) fa + ^ / / alü<L·'0(u),u')fa mm/ \u,u)',rn,rn , (^})fa 

tmm > JO imm> JO JO 

where \(j,m,m'(<f))fa — symb\bj,m,m')fa and \u>,uj' ,m,m' ((/))) fa = symb\uj,Lü' ,m,m!)fa. 

c- The quantum states p are measured by the observables just defined, leading to the mean vàlues of these 
observables; in the usual notation: (O)p = TrffiO). We can conceive that mean vàlues as the more primitive objeets 
of the quantum theory (see [52]). These mean vàlues, generalized as in paper [24] and symbolized as (p\0), can be 
considered as the result of the action of the linear functionals p on the observables of the vector space O. Then, 
p G ScO , where S is a convenient (i.e. satisfying eqs.(45) and (46) below) convex set contained in O , the space of 
linear functionals over O. The basis of O' (that is, the co-basis of O in each D fa) is {(uj, mm'\fa , (uu', ram! \ fa'], and 
it is defined in terms of its functionals by the equations 

(uj,m,m'\fa\r],n,n')fa = 5(u> - ■n)5 mn d m > n >5 i j 

(u!,u)',m,m'\fa\r],r]',n,n , )fa = S(w - rf)5{uj' - r)')5 mn 8 m > n >5ij (42) 

and all other (.|.) are zero. The orthogonality in ... is a consequence of eqs.(37) and (41). Lct us observe that 
(u>,u)', m, m'\fa = \uj,m)fa(oj',m'\fa but (u),m,m'\fa ^ \oj,m)fa{uj,m'\fa. 18 Then, a genèric quantum state reads 



17 The component Òs is called singular because it contains a hidden distribution S(lo — lo'). In fact, it can be obtained from 
the regular part by making 0(u,u% iTnm > = 0{uj)fa mm ,5{ij0 - uj'). 
18 If (ai,m, m'\fa — \üj,m)fa(üj,m'\fa, it is easy to show that a divergence appears. 
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p OO p OO p oc 

P=J2 dwp{w) (j>imml {u,mm l \ < t >i + / dLJ ^V(w,w')^ w (^',mm'| 0t (43) 
,Jo ,Ja Jo 

imm imm 

This is also the case for the corresponding classical analogue p(4>) of p, 



JO Jo Jo 



(44) 



Each of the terms of the sum J2i can be considered as a tcrm of a decomposition, where each of the p^ t (</>, t) = 
symbp^^t) does not vanish in the corresponding domain D^ >i C X^. 19 As before, the first term in the r.h.s. of 
eq.(43) will be called ps, the singular component, and the second term will be called pn, the regular component. 
Functions p(ui, v')<p imm ' are regular (see [24] for details). 

Going back to the decomposition (44), we impose the following conditions. We require that p = p, i.c. 



P( w > u') é , = p{J, u)^ m , m (45) 



and that p(w)^. mm , be real and non-negative, satisfying the total probability condition 

/•OO 

(p\T) = £ / dupiw)* = 1 (46) 

im J 

where I = J„°° du) ^2 im \ui, rn) ( p i (ui, ra\ ( p i is the identity operator (24) in O represented in each D^. Eq.(46) is the 
extension to state functionals of the usual condition Trp 1 = 1, when p is a density operator. Thus, from now on, 
Trp = (p\I). For these reasons, p belongs to the already defincd convex set <S C O' . The time evolution of the 
quantum state p reads 

/•OO /'OO /'OC 

p(t)=J2 diüp-^) (üj,mm'\^+ / ^ ^'pM^ m( e i(uV)t/í (^>m'| fc (47) 
/Jo / Jo Jo 

imm' imm' 

As we have already said, at the statistical quantum level we can only measure mean vàlues of observables in quantum 
states 

/>oo /•OO rOO 

(Ò) m = (p(t)\0) = + E / <W dLü'rt^ ^-"'^O^M^mm' 

, JO t jo jo 

imm' imm' 

(48) 

From eq.(16), the classical analogue has exactly the same form 

/'OC 

{0{<t>)) p(M = {p{^,t)\o{4>)) = J2 / <^p^)^ mm ,o{Lü)^ mm , + 

, JO 

imm 

pOC pOO 

+ J2 du; du/^^^/e^-^^O^u/kw (49) 
, ,Jo Jo 



Both decompositions are vàlid in each D ( j >i . If we take into account that 0{oj,ijj') and p(w,w')^ imm / are regular (as 
regular as needed to use the Riemann-Lebesgue theorem, i. e. 0(u>,u>') p(u> , uj') ^. mm i eLifw- ui'), see [24]), we can 



19 Considering for a moment the larger domains T)^, of course there are other charts defined in other domains and, 
therefore, other D'^ k C 2^, fc - But since p(4>,t) = symbp(t) is defined in the whole M, the T>^ >i , V'^ k are just local charts for 
which the same function p(<j>, t) = symbp(t) is globally defined in phase space. Moreover, at 4> £ T>$ i n V^, , any pair of charts 
can be C T — smoothly connected in the sense that all their elements can be smoothly connected among each other. The same 
argument can be applied to the partial decompositions of pn{<j),t) = symbpn(t), ps{4>,t) — symbps(t) (i.e. the first and second 
terms of the r.h.s. of the last equation), since both regular and singular parts are also globally defined in phase space. 
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take the limit t — ► oo and use the Riemann-Lebesgue theorem. As the result, we see that the fluctuating-regular part 
vanishes and we arrive to the weak (quantum and classical) limits 

W lim p(t) = = % = V / cL;p(w,p), m ,(w,m,m'|^ 
t^°° , Jo 

imm' 



W lim p(0,í) = p s (0) = P*(</>) = V / dujp{(j,p), m ,(w,m,m',((/>)|^ (50) 

, Jo 

imm' 

where (o>, m, m', (<^)|^ 4 = symb(uj, mm'\ ( p i and = symbp* are defined in A'l and the integral in eq.(50) is decom- 

posed in diffcrent ways at each D t j >i . Since only the singular diagonal terms remain, we have obtained decoherence in 
the energy variable lo. Precisely, any quantum state weakly tends to a linear combination of the energy diagonal states 
(ui, m, m'\ ( f >i (the energy 'off-diagonal' states (u),u)', m, m'\fa are not present in /?»). This is the case if we observe 
and measure the system evolution with any possible observable of space O. Therefore, from the observational point 
of view, we have decoherence of the energy levels in spite of the fact that, from the strong limit point of view, the 
off-diagonal terms never vanish: they just oscillate since we cannot directly use the Riemann-Lebesgue theorem in 
the operator equation (47). 
Important remarks 

i. - It may be supposed that decoherence takes place without a coarse-graining. It is no so: the choice of the àlgebra 
A among all possible algebras (see under eq.(39)) and the systematic use of mean vàlues (0)^ t ) — (p(t)\0) (eq.(48)), 
restrict the availablc information and produce the effect of a coarse-graining. In fact, we can define the projector 
II = |0)(po|, w ith \Ò) £ A and (po\0) — 1, that projects (p(t)\ as (p(t)\U — {O)p^(p \, and translates everything 
in projectors language: we obtain, from eq.(50), linií^oo^í) II = (/5*|IT. This projection will obviously break the 
unitarity of the primitive evolution. In this way we could develop a formalism closer to the usual one. See a detailed 
cxplanation in [13] and [34]. 

ii. - Theoretically, decoherence takes place at t — > oo. But, in practice, decoherence appears at a decoherence time, 
as we have defined in [53] : the decoherence time can be easily computed from the poles of the resolvent or the initial 
conditions density in the complex extension of the H spectrum. Trivial H (e.g. free particle H) and trivial initial 
conditions (e.g. zero temperature ones) do not have poles and the decoherence time is infinite. This means that, to 
reach equilibrium in a finite characteristic time, H must be non-trivial (e.g. the sum of a free Hamiltonian plus an 
interaction Hamiltonian) and/or the initial conditions must be non-trivial (e.g. T ^ 0). For details, see [54], where 
decoherence times are estimatcd in 10~ 13 — 10~ 15 s for microscopic bodies, and 10~ 37 — 10~ 39 s for macroscopic bodies; 
for a thermal bath our rcsults coincide with those obtained by the einselection approach. 

C. Decoherence in the remaining variables 

Having obtained decoherence in the energy levels, we must consider decoherence in the other dynamical variables 
( j )i i of the set of local CSCOs we are using. We will call these variables 'momentum variables'. Since the expression of 
given in eq.(50), only involves the time independent components of p(t), it is impossible that a further decoherence 
process eliminates the off-diagonal terms in the remaining N dynamical momentum variables. Therefore, the only 
alternative is to find the basis where these off-diagonal components p{io)^ im m' vanish at any time. 

Let us consider the following unitary change of basis 

\u,p)4>í = y^ J U(u) mp \uj,m)j, i (51) 

m 

where p and m are shorthand notations for p = {p\, ...,pn} and m = {m\, ...,tojv}, and [[/(w) -1 ] = U(u) pm . We 
choose the new basis {Jc^p)^} such that it verifies the generalized orthogonality condition 

(w,pU|w',p% 4 = 5(w -uj')5 pp , 

Since p(w)^. = p(^)(j>i, h is possible to choose U(ui) in such a way that the off-diagonal parts of p(w)^ iP p' vanish, i.e. 

p(u)4>ípp' = P<t>i(u) P S pp' ( 52 ) 
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This means that there is a final local pointer basis in D í f >i for the observables, given by {\ui,p 1 p l ) ( j >i , \üj,u)' jPjp')^} 
and defined as in eq.(41) but now with the p. The corresponding final pointer basis for the states, {(^,P,p'\^i, 
(Wjü/jPjP'I^}, diagonalizes the time independent part of p(t) and, therefore, diagonalizcs the final state 
Now, we have diagonalized the p{u)^. mm , in m and m', obtaining 

roc 

W lim p(t) = as = p* = V / dujp^iuj) J(J,p,p\^ 



_ roc 

W lim p((j>,t) = p s {4>) = p*{4>) = Y] / düjp^iuj) (u,p,p, (çí>)| 0i (53) 

ip " 

Here we are using a local pointer N + 1— CSCO {H, P^i, P^ív} at each D^, where the P^/ are 

p<t>ii=y~\l dwy^p^/Miw,p,p)^ (54) 

i •'o p 

and their classical anàlogues 

/>oo 

p {<t>)<t>ii = y] / ^y"í , 1 /(^)l·,p>p(</ , ))0 1 (55) 

where la;,^,^)^ = Pl<^>i or simply {jw,;}}^} is the foca/ pointer basis in D^; so, we can write eq.(40) in 

this new basis (see eq.(58) below). 20 Now all the operators and matrices involved are diagonal, and decoherence is 
complete. We can define all the observables O of eq.(40) in this new local pointer basis. 

Since in the limit H — > we usually have P with continuous spectra, instead of the last equations we would have 
the natural anàlogues of eqs.(53) (see [31] and [33] for details) 



W 

t 



lim p(t) = as = P* = V / du ( dp N p(uj) 4>z (uj,p,p\ c f )t 

i JO JpcD+ { 



W lim p{4>,t) = p s {4>) = P*{<j>) = V / du dp N p(uj) 4>z (uj,p,p,((f))\ c j >z (56) 

'-* 00 i Jo JpcD+ t 

In the next section we will consider the classical limit and, then, we will only use continuous spectra and equations 
like the last two: 21 so we will re-write some equations in the new basis for the sake of completeness. 



20 The complexity of these formulae demonstrates why it was so difficult to define the pointer basis in a general case. As we 
can see, the pointer basis depends on H and the initial conditions, but there are some cases (see section IV) where it only 
depends on H . 

21 If we use the Heisenberg picture, the A would become diagonal. So, heuristically 

lim (p,\[À(t),B}) = lim Tr(p,À(t)B - p,BÀ(t)) = 



Tr(p*À*B - p t .BÀ«) = Tr{p r À*B - A,p,B) = 

where A, is the diagonal weak limit of A and, therefore, commutes with the diagonal p*. As a consequence, the evolution is 
(heuristically) weakly asymptotically abelian ( [51], Def. 4.11) since, in the limit t — > oo, A can be considered commutative. 
Therefore, a quantum system with continuous spectrum is weakly asymptotically abelian. 
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IV. THE CLASSICAL STATISTICAL LÍMIT 



A. Quantum and classical operators 

a.- From now on, we will consider a system from the point of view of the local pointer completo set of N + 
1— commuting observables {H^P^i, P^n}, defined in eqs.(35) and (54). As above, to simplify the notation wc 
will just call {HjPfc} the set {H,P^ t i, P^n}- Thus, we will consider the orthonormal eigenbasis {\w , p) fc} of 
{HjPfc}, 22 and write the Hamiltonian and P as 

H = / dp N / wl^p^^pl^dw Pfa = I dp N / p^.p)^^^^^ (57) 

i JpeDj,. JO JpeD^ Jo 

Furthermore, we will consider the àlgebra A of the operators (40), which now read 



= V/ dp N O 0i (w,p)|w,p) 0i (w,p| 0í GÍu; 

„• JpeD^. Jo 



^ / / dp N dp' N 

IpeDj,. JptD 4 , t JO JO 



OO />00 



^ I i dp N dp' N I I O rt>i {LO,J,p,p')\^,p) 4>i {LO,p l \^(L·(L· l (58) 



As before, the first term in the r.h.s. will be called Òs , the singular component, and the second term will be called 
Or, the regular component. Also as before, functions 0^(0), u>' ,p,p') are regular (see [24] for details), [H,Os\ = 0, 
Os € Cs, where Cs is the singular space, Or e Cr, where Cr is the regular space, and A — Cs®Cr. The observables 
are the self-adjoint operators of A, and they belong to a space O. 

b - Let us now consider the Wigner transformation of these objeets. The operators of Cr are regular; so, their 
transformation is obtained as explained in Section II. Thcn, wc have to consider only the singular space Cs, the space 
of the operators that commute with H. This is not a regular space of operators on a Hilbert space TL as Cr, since 
it contains a hidden 5(lü — u>' ) (see eq.(39)), but the mapping symb given by eq.(4) can also be well defined for the 
observables in Cs- In fact, from eq.(58) we know that 

Os =51/ dpN ^ 4 (w,p)|w,p)^(w,p|^dw (59) 

i JpíD 4>% Jo 

If we consider, as usual, first ( j >i as a polynomial, and then < f >i as a function of a certain space where the polynomials 
are dense, 23 by using eqs.(57) we can conclude that 

Os = ÇO*(íï,í£) = Çds* (60) 

i í 

But, when f,g commute, as the members of the CSCO do (see eq.(12)), 

symb{f,g) = (/*<?)(<£) = f(4>)g(4>) + 0(ft 2 ) (61) 
Then, by means of the same procedure as before and eq.(7)), 

symbÒs = O s {4>) = Ç °4>i( H (<t>)> P 4>M) + = Ç s V mb Os^ (62) 

i i 

where -ff (</>), P ( j >i (<j>) can be computed as usually (sec [31] for details). In this way, we have succeeded in computing all 
the symb of the observables of Cs up to 0(h 2 ), which are just the 0^ i {H{<j>), P(4>)), and we have defined the mapping 



22 We are using the 'final pointer basis' of section III. C. Below we will write all the formulae in this basis. 
23 These polynomials have several variables, but there is no problem since all these variables commute. 
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symb : C s -> £s 9 symbO s = O s (0) = ^ {H(<t>), P^ (0)) + 0(h 2 ) 



(63) 



Moreover, since decompositions or are arbitrary (because they depcnd on the initial conditions of Section 
III. A), from eqs.(60) and (62) we obtain (up to 0(h 2 )) 

Ò s ^=0^(H,P^), OstM = symbd s <t> i =0 4>i {H{<l>),P<i, i {<i>)) (64) 

Lct us observe that, if 0^(10, p) = 5(u> — oj')5(p — p'), we have (also up to 0(fi 2 )) 

symb^^p')^'^'^ = S(H(^)-u;')S(P^M-p) (65) 

an equation that we will use below. 

Summing up, from eqs.(2) and (63) we have defined a classical space A q — Cr © Cs and a mapping 

symb : A q symbO — 0{<f>) (66) 

where eqs.(9) and (10) are also vàlid. Thcn, we can repeat what we have said below eq.(10), but now for the àlgebra 
A q defined as in this section, with its regular and singular parts. 

If now we take the limit fi — > 0, we obtain A q — > A, where A is the usual àlgebra of observables on phase space. 
Thcn, in this limit we have a correspondence A.—* A. However, even if this limit is well defined and can be considered 
as the classical limit of the àlgebra of operators, it is only the limit of the equations of the system, sinec these are a 
consequence of the àlgebra. Therefore, this is just a 'formal' limit. The limit fi — > will be completely studied when 
we deal with the state space. 

For the sake of simplicity, from now on we will systematically climinate all the 0(fi 2 ) from the equations and call 
the A q just A. This is a rigorous simplification. In fact, when fi = wc can makc the 0(fi) = everywhere since, 
from cq.(7), when fi = we have expO = 1 in that equation; in other words, the lim^^ is continuous. 



B. Quantum and classical states 

a.- Let us remember that \w,p) cí , i = |w,p,p) 0ï = \u:,p) 4>i (uj,p\ 4>i and \uj,lj' ,p,p') lt>i = as in cq.(41). 

{\u),p,p')fa} is the basis of Cs and {\ui, o/,f>,p')<£«} is the basis of Cr. Then, cq.(58) rcads 

d = T, d = e / d P N r o^p^p^du 

i J PíD^ Jo 

p p pOO pOO 

+ E / / dp N dp' N / / 04 t fau , ,p,p')\u,u>',p,pf)4 i <L·<L·' (67) 

i JptD 4 , t JpeD 4>l· JO Jo 

Since the states are functionals over the space A — Cs © Cr, let us consider the dual space A' = C' s © C r. We will 

call {{^,p\4>i} the local bases of C' s and {(w, w',p,p'\fa} the local bases of C' R . Let us remember that (oj,u)',p,p'\ ( j, i — 
\uj,p) íjlz {uj' iP 1 ^ but (w,p|^ 4 ^ \oj,p) ljlz {uj,p\ < p i . Moreover, as in eq.(42), 

{uMïiW ,p')<l>3 = s (u - u')8 N (p - p')5ij (u,a,p,s\^\u}' \a' ',p' >')^ = 5 ( UJ ~ u')5(ct - a')5 N (p - p')S N (s - s')5 tJ 

(w,o-, UJwV.í/.s')^ = (w,cr,p,s|0 1 |w',cr')0 i = (68) 
Then, a genèric functional of A' reads 

P = EÍ dpN í / 9 <fc( w >P>)( w >ï , k dw 

i JpnD 4 , i Jo 

+ E / / dp N dp lN ^ r P0I (w, lü', p, p') (üj, w',p, p' I fràudu' (69) 

i JpeD^,. JpeD 4 , z Jo Jo 
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Likc functions ( p i (lo, u>' ,p,p'), functions (w, u>',p,p') are regular and have all the mathematical properties necessary 
to make the formalism successful (see [24]). Moreover, the p must be self-adjoint, and their diagonal p^^.p) must 
represent probabilities; thus, J °° p c / >i (uj,p)duj — 1 (as in eq.(46)) and, most important, 

P4u (w,p)>0 (70) 

The p with such properties belong to a convex set S, the set of states. Also, as in eq.(48), 

P P oc 

(p\Ò)=y2 / P^ i (^,P,)04> i {uJ,p)(L·dp N 

i JpeD 4 , z JO 



P P P OO POC 

+ V / / / / P4>t ( W , U', P , P')0^ (W, Uj',p, P ')dLüdLü'd N p d N p' (71) 

• JpeDj,. JpeD^ Jo JO 

b.- Since Cr and C' R are spaces of operators on a Hilbcrt spacc H, the symbol for any p R £ C R is defined as in 
eq.(15). 24 From this definition, eq.(16) can be proved for the regular parts (see the demonstration in [43], eq.(2.13)): 

(pr\Or) = (symbp R \symbOR) = V / d^ N+1 ^ p^ R ^)0^ R ^) (72) 

i Jd 4-í 

Then, in C R and C' R all the equations are the usual ones (i.e. those of papers [43] and [44]). 

Let us now consider the singular dual space C' s , the case not treated in the bibliography. In this space we will 
define symbps as the function on M that satisfies an equation similar to eqs.(16) or (72) for any Os S Cs, namely, 

(ps\O s ) = (symbps\symbO s ) 

precisely, 

V/ r P(l>i (uj,p,)0^{uj,p)dwdp N = V / d^W+Vp^sMO^sW (73) 

where the unknown density function ps {(fi) — symbps can be decomposed as 

symbps = Ps(4>) = ^2 PfasW) ( 74 ) 

i 

in each -D^. Thus, since we know p^,i(uj,p,), 0^(0), p), and O^siffi), we can compute Pfcsi'fi) to obtain ps{4>) = 
symbps- Now (as suggested by eq.(62)), ps, being time invariant, must be a function of the constants of the motion; 
therefore (as in Subsection A) its Weyl-transformed ps ((fi) must be endowed with the same property, but now in the 
classical case. Since the {H((fi), Pfcicfi)} are locally a complete set of constants of the motion in involution, we must 
have 

Prt>i (cfi)=F(H(<fi),P4>M) (75) 

We will find the function F. The system has a local pointer CSCO of N + 1 operators and the dimension of its 
phase space is 2(N + 1), i.e. it is locally an integrable system. 25 Then, we can define locally at D ( / >i the action angle 
variables (6°, 6 1 , N , J°. , J^., J^.), where J°. , J^., would be just H, Pfai, P^n (multiplied by adequate 
constants in such a way that the 6 l P. variables belong to an interval < 9 l P < 2ir in the integrable case). We will call 



24 We repeat that, in the case of states, we must add a new factor (2nh) ( JV+1 ) to definition (4) in order to maintain the usual 
normalization of p((j>). 

25 We have discussed this fact in detail at the beginning of section III. The constants J are global or isolating in the case of an 
integrable system, but not in the non-integrable case. Nevertheless, they are locally defined. Moreover, we will only consider 
the cases where action-angle variables can be locally defined. 
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'J' just the , H } P^'. Thus, we can make the canonical transformation <fi a — > ...,9^.,H, P^i, ■■■,P^ >í n, an d we 

obtain 

d0 2 ( N+1 ) - dg^dp^ - d0% +1) dHdP£ (76) 

because the Jacobian of a canonical transformation is ±1. Since the transformation must keep the mètric of eq.(3), 
we can take 1 with no loss of generality. 

In order to compute the l.h.s. of eq.(73), we must know how to integrate functions f(H,Pfa) = f(H, P^, ...,P r j >i ) 
which are just functions of the constant of motion, precisely, 

/ d4> 2 < N+ »f{H,P <l·t )= í de^ +1) dHdP^f(H,P^) 



= í dHdPgC^{H,P^)f[H,P^) (77) 

where we have integrated the angular variables ■»,0$ i and obtained the configuration volume C^, t {H,P^ ( ) of 

the portion of the hypersurface defined by (H — const.,P c / )i = const.) and contained in D^ t . So, from eqs.(73) and 
(77) we have that 



peD^ JO 



p<jn(u,p, (Lú,p)divdp 



N 



= j dHdP^C^{H,P^)p^s{H,P^)0^s{H,P^) (78) 

for all 0^(11, P 0i ) = S4>i (H, P 0i ) (see eq.(64)). The last equation defines p s<t>i {H,P) = ^-p 0i (ií,P) for e X> 01 , 26 

but not for </> e AlvX^ 27 ; then, as in the case of Os^ (</>), we will consider that ps<j>i{4>) = for cj> G MXD^ and that 
they are defined all over M. (this causes no problem because Os<t>i{<)>) is multiplied by ps<t>i(<f>), an d Osfcicf)) has this 
property). In this way, we can arrive from eq.(74) to our final result 

Ps{4>) = PM = E r À P s POi (H(<P), P* t m (79) 

Now, from eq.(70) we obtain that 



26 In the integrable case, where there is just one p(H,P), it would be p^^HjP) — P) an d tne results of paper 

[26] would be reobtained. In fact, integrating over a torus in the 8 we have (2iv) N+1 p(H, P) — JT C c j >i (H, P)p^ i (H, P). 

An example to fix the ideas: let us consider the harmònic oscillator and the plane q, p in radial coordinates 9, H . Let us define 
two D4>i\ Di with < 6 < Q(H) and D 2 with Q(H) <0<2n, where Q(H) is an arbitrary function. Then, 

p{4>) = piWiW + /»(4>)W4>) 
If p(0) = p(H), by integrating over the 6 we obtain 

fB(H) r'2n 

2np(H) = pi(H)h(<t>)d0 + p 2 (H)h(<j})de 

Jo Je(H) 



= Pl (H)Q(H) + (2tt - Q(H))p 2 (H) = 



P i(H)Ci(<l>) + p2(H)C 2 (4>) 



namely, the equation (2ty) n+1 p(H, P) — JT C ( p i (H, P)p<j> i {H, P) for this particular case with dimension N + 1 = 1 
27 We will forget the joining zones J 7 ^ and F ( p i since we have already proved that, when S — > oo, their influence is irrelevant. 
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ps(4>)=p.(4>)>0 



(80) 



This means that the Wigner transformation of the singular part can be considered a density function since it is 
non-negatively defined (of course, this is not the case for the regular part). 

Always working in the domain V t j >i and making p^^uj.p) = S(u) — ui')S N (p — p'), we also have 28 

symb{J,p', {<f>)\^ = ç^p^ S {H{4>) 6W (P{4>) ^) (81) 

c- From eqs.(56) and (81) we obtain 

Ps{4>) = P*{4>) = E / dp í°° P4>i (uj,p) - l 5(H^)-u)5^{P{cj > )-p^)tL· (82) 

i JptD^ Jo L· <M- H ' / >J 

The continuity of the function p*(<fi), when it goes from one D ( f >i to another D ( j >j (i ^ j), is demonstrated in Section 
V. Therefore, we have obtained a decomposition of p*{4>) — ps(<t>) in classical hypersurfaces (H = uj, Pfafy) = 
Pfc), containing classical trajectòries, summed with different positive weight coefficients p<f >i (u),p)/C < i >i (H,P ( f >i ), and 
represented in different ways in each domain D^, but still with the same interpretation as in the integrable case. 29 

d.- Since now we know how to deal with the singular part, we have defined the mapping of the quantum space of 
statcs A' on the 'classical' space of states A' 

symb : A> -» A' (83) 

In the limit H — > 0, eqs.(9) and (10) are always vàlid; then, it might be supposcd that we have arrived to the classical 
limit for the states. But this is not so because, in general, even for S very big (or H very small) the obtained p(<f>) 
does not satisfy the condition (see also the Appendix of rcf. [31]) 

p{4>) > (84) 

This is duc to the fact that the regular part is still present and this part does not satisfy the last condition (on the 
contrary, from eq.(80) we can see that the singular part satisfies the last inequality). As a consequence, p(<j>) is not 
a density function and, therefore, the mapping (83) is not a mapping of quantum mechanics on classical statistical 
mechanics. This mapping does not give us the classical world, but a deformed classical world where 'density functions' 
can be negative. In other words, when h — > the isomorphism (83) is a mapping of quantum mechanics on a certain 
quantum mechanics 'alia clàssica', namely, only formulated in phase space M but not satisfying cq. (84). This 
clearly shows that ^ — > is not the classical limit. In order to obtain this limit, we have to introduce decoherence, as 
previously studied, both at the quantum and the classical level. 



28 In the chaotic, homogeneous, ergodic case, we have a N + 1— CSCO with just H and, classically, just H as a constant of 
motion. In this case (see [55], p.247), 

S(H(<j>)-E) 



f°° 1 
ps{4>) = P*{4>) = j o S(u — E)-——S(H((f>) — E) 

29 We can verify the normalization: 



fdqdp5(H(cj>) - E) 



j Ps {4>)d^ N+1) = Ç/ ps^)dHdP»d8^ = 
Ç j dHdP£ Ps (H,P^) = 1 
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C. Time evolution and decoherence 

As we have seen, the only thing that prevents us from having a good isomorphism (83) is that the regular parts 
do not satisfy condition (84). But we know from eqs.(50) or (56 ) that, for t — > oo, the regular part vanishes and 
only the singular part remains, which does satisfy this condition. As a conscquence, after decoherence and fi — > 
(that is, the elimination of all the 0(fi 2 ) that we have omitted), we finally obtain the classical statistical limit since 
the classical densities obtained obey all the laws of classical statistical mechanics. In fact, as we will see in the next 
section in detail, eq.(82) shows that these distributions are the result of classical point-like-states moving in phasc 
space and following classical trajectòries. The usual classical limit is obtained by choosing one of these trajectòries; 
we will explain this procedure in the next section. 

V. THE CLASSICAL LIMIT 

From what we have learnt above, we can explain with more detail the three steps involved in the classical limit, 
presented in the introduction and shown in the following diagram: 

Quantum Mechanics — (decoherence) — ► Boolean Quantum Mechanics — (symb and fi — > ) — > 

Classical Statistical Mechanics — (choice of a trajectory) — > Classical Mechanics 

Let us comment these three steps: 

i.- Quantum Mechanics — (decoherence) —> Boolean Quantum Mechanics. Decoherence transforms non-Boolean 
quantum mechanics into Boolean quantum mechanics 30 since it eliminates the off-diagonal terms, as we have shown 
in eq.(56). 

'ú.-Boolean Quantum Mechanics — (symb and fi — > ) — > Classical Statistical Mechanics. The Wigner trans- 
formation symb and the limit fi — ► are defined with no problems in the singular part remaining aftcr decoherence. 
They map Boolean quantum mechanics onto classical statistical mechanics: this is what we have essentially shown 
above. Our demonstration culminates in Section IV. C, where we have proved that the transformed quantum Boolean 
states are really positively defined densities. From eq.(82) we also know that these densities are the sums of densities 
strongly peaked on the classical hypersurfaces defined by the constants of the motion H(<j>) = u, P^i{<f>) — Pfc- In 
the next step we will see that such classical hypersurfaces contain classical trajectòries averaged by the coefficients 
P0i(w,p). 

iii. -Classical Statistical Mechanics — (choice of a trajectory) — > Classical Mechanics (Localization 
or Actualization). After step (ii), we are still in classical statistical mechanics but not in proper classical me- 
chanics. To perform the last step we have to pass from classical densities to classical trajectòries (i.e. to consider the 
localization effect 31 ). For this purpose let us observe that, after the two first steps, the formalism of Boolean quantum 
mechanics is isomorphic with the formalism of statistical classical mechanics: 

• For the observables: After symb and fi — > 0, we obtain the correspondence symb : A ~ A (see Section IV A), 
namely, 

A^H^) ~ A^H^P^)) 

• For the states: After decoherence, symb and fi — > 0, again symb : A' ~ A' (see Section IV. B), namely, 

PtMhò "P^mt), (</»)) >o 



Namely, quantum mechanics in the local CSCO {H, P} using only diagonal states. 
See [47], Chap. 4, for a different view. 
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and the states p*(H,P) and p^(H((f>), P(<j>)) are time invariant: 



p*{H{<t>),P{<t>))=Y J fduí d P N ^^6(H(iï-co)5 N (P <P M-P<p t ) (85) 

Moreover, since AAAB > §|(L4, B]) p \, in the limit fi — > there are no uncertainty relations and the algebras A and 
A can be considered commutative (remember that, according to the uncertainty principle, fi — ► has the same effect 
that [A, B] =0). In other words, in the limit fi — > all the picture is classical in such a way that the trajectòries, 
containcd in the hypersurfaces H{4>) = ui, P^cf)) — p^, could be interpreted as real classical trajectòries. However, 
the S(H((f>) ~ ui)S N (P0í(</>) ~ P<t>i) still represent states strongly peaked around these hypersurfaces (but not around 
trajectòries). Therefore, if we want to obtain an equation like (85) but clearly showing the classical trajectòries, we 
have to introduce the initial conditions of each trajectory. 

Let us consider a classical trajectory in phase space Ai = M.2(n+i)i expressed in the classical coordinates 
(r, 0^, H, Pfc), where r is the coordinatc canonically conjugated to H and the 9 ( p i are the coordinates canonically 
conjugated to the P^. The constants of the motion in involution are {H, Pfc}; but, for conciseness and generality, 
let us consider that the constants of the motion in involution are {ü^} with conjugated coordinates {A^}, and that 
H = HiTlfa). From the von Neumann-Liouville equation in the Heisenberg representation, 

,r,f = { A,H\ 



we obtain 

dA{<t>) 
dt 

Then, the Hamiltonian cquations in the limit fi — > read 32 



{A,H} mb = {A,H} pb + 0(ti 



2\ 



dA+ t dH dU^ 8H 

^r = ^- = ^ (n ^ ) = consí · ; — = -dAjr° (86) 

The classical trajectòries are 

a^ (t) = Af t + n 0( (n 0i )t, n 0í = njj = const. (87) 

where the A^ and are integration constants. A distribution strongly peaked on this trajectory reads 

and a general classical distribution evolving according to the motion (87) rcads 33 

Pc{tA) = E / pf^A^^SiA^t) - - 0^(n^)í]x 



SiTL^-TL^d^A^d^TL^ (88) 

where p£' (A^ , ) is a genèric classical coemcicnt (undcfincd up to now ) . If we want that this density (evolving 
according to a Frobenius-Perron evolution, sce [56]) be an equilibrium density, we have to eliminate the variable t. 



32 These equations correspond to the system of differential equations (3.1) of [56]. 

33 If the evolution S f of [56] were the (87), the corresponding density would be f{t,x) = P t f(x) (see [56], eq.(3.2)) where P f 
would represent a Frobenius-Perron evolution. Moreover, it is easy to show that p(t, 4>) satisfies the Liouville equation. 
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For this purpose, it is suficient to choosc the initial distribution p^\A^ ,11^) as just a function of ü^, namely, 
p^(IT^), which is still a frec function that wc can usc to rcprcscnt diffcrcnt pc (</>)■ Thcn, wc obtain 

pc(^) = £ / pWtfgVV* ~ n^^nw (89) 

since, for any fixed t, we have 

E / - A S? ^ ^.(n.Jí)^ 1 ^ = i 

Going back to our primitive variables, eq.(88) reads 

Pc{4>) = E/ P^^,p)S{H(4>)-w)S{P 4ti - P4>i ) 



xS(r(<f)) - r - wí)í(^ 4 (0) - 6»^o - p<p z t)dLjd N p^drodO^o (90) 

while eq.(89) reads 

PcW = E / pfffaPMHtt) - ^(^ -P^d^p^ 
Thcn, from cq.(85) and making the undcfincd coefficient p^\uj,p) — ç*^"^ ) i we have 

pM = pc{4>) (91) 

The function pc(4>) can be interpreted as the equilibrium density of a Frobenius-Perron evolution of particle-like states 
(r, 0,^, ií, -P^), as i/these states would move in phase space Ai — M-2(n+i) according to the classical motions (87). 

However, each term of the sum J2i 01 eq.(90) is vàlid in the chart V t j li (Dfc C T)^). In a diffcrent chart T>^ (D^ 
C T>j,j), the equation is also vàlid and, then, at <j> <E T> ( p i n we have 



xS(t((I>) - t - ut)S(9 < f >i (4)) - 6>0 iO -pfatfdwd p^drodA^o = 

x5(t(4>) -t - Lüt)5(6^((f>) - O^o - p r f, j t)dLjd N p (j , j dT dA (j , j0 

Here it is worth emphasizing that the trajectòries ií = w, P^^cj)) — p^, r((f>) = to + wí, 0^(4)) — 9^0 +P4>it in chart 
2?0 i are continuously connected with those H — uj, P ( p j (<f>) = p r j >j , t(0) = to + u>t, 9 r j >j (<fi) = 9^0 + p^t in chart T> r f >j , 
because these charts are not genèric but constructed using the solution of eqs.(19), (23) or (86). Since D c / >i C T> ( p i 
and D t j )j C T> < j >j , the same holds for the trajectòries going from D ( f >i to D < f >j . Thus, the continuous connection follows 
from the fact that one and only one solution of the trajectory equation passes for each point of M. (and, therefore, 
for each e n V^). 

Summing up, we have obtained a decomposition of p*{<f>) = ps{<fi) m classical trajectòries H = u, Pfa^tf) = Pfa, 
t{4>) = t + ut, 8 < p j (<fi) = 6*0 jO + P<f>j t, summed with different weight coefficients p^ (u>, p) / C l p i (H, P t f >i ) and represented 
in different ways in each domain T> l p i , but still with the same interpretation as in the integrable case. Moreover, as 
announced in Section III. A.c, we see that chart T> c f >i is continuously connected with chart "Dfa, for any V^, "D^. 
Therefore, we have finally obtained the classical limit to the extent that we have described each one of the classical 
trajectòries. But, since from the very beginning our system was a non-intcgrable one, we have obtained the classical 
limit of a non-integrable system, where the tori are broken and the trajectòries are potentially chaotic trajectòries. 

Finally, we must remark that: 
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• Each one of the described processes, decoherence, route to macroscopicity, i.e. h — > (e.g. the macroscopicity 
obtained when the two rays of an Stern-Gerlach experiment gradually separate), and eventually localization (e.g. 
by a localizing potential, see [26], Appendix A), has its own characteristic time; in particular, the decoherence 
time is computed in [54]. 

• We have explained the classical limit as if each process (decoherence, macroscopicity, and localization) took place 
one after the other, only for didactical reasons. But this is an oversimplified picture of the phenomenon. In fact, 
this may be not the case if the different processes overlap. Considering that they have different characteristic 
times, there are different possibilities according to the order in which the processes finish. 

VI. PARTIALLY NON-INTEGRABLE SYSTEMS: THE TRACE OF THE NON-INTEGRABLE PART 

In this section we will consider some well-known definitions, find some relations with the concepts of statistical 
physics, and discuss the trace of the non- integrable part of the system. All the definitions are formulated in the classical 
language, but they become quantum definitions if we translate them with the Weyl-Wigner-Moyal isomorphism. In 
some sense, that we will precise elsewhere, these quantum definitions correspond to a quantum chaotic hierarchy. 

The classical constants of the motion H(<j>), P\(4>), Pn{4>) can be rigorously classified as [57], [58], [59]: 

i. - Global or isolating constants of the motion, that we will call 'H\ 

H (<f>) = H(<f>), ffi(0) - Pi (</>),..., Ha (</>) - Pa(<I>) (92) 

The constants of the motion are global when the conditions 

H a = Pa , a = 0,...,A (93) 

dcfine, for each set of constants (po, ...,pa), a global sub-manifold A4(po, ...,pa) (a torus in the bounded case) of phasc 
space where the trajectòries are confined. 34 The dimension of M. (po, ■■■,Pa) is 2(N + 1) — (A + 1) = 2N —A + 1. 

ii. - Local or non-isolating constants of the motion, that we will call 'J^', 

JfriW) = P^a+i(4>), J^ t N-A(<f>) = P*M<t>) ( 94 ) 
The constants of the motion are local when the conditions 

J<f> if 3 = Pfrp+A 0=\,..,N-A (95) 

do not define any global sub-manifold, but define just a local sub-manifold at D ( p i . 35 The J^p are local momentum 
coordinates. 

When A = N, we say that the system is integrable; when A = we say that the system is non-integrable (even if 
we will always consider H globally dcfincd), if A < N we will say that the system is partially non-integrable. Let us 
consider these three cases in detail. 



A. Integrable systems 

In this case, all the P are isolating constants of the motion H, and there are no J ( f >i . Then, condition (93) foliates 
the phase space with submanifolds M(po, ...,pn) of dimension N + 1, labelled by the constants (po, ...,pn)- These 
submanifolds are tori when the system is endowed with action-angle variables (e.g. when phase space is bounded). 
On these tori, the motion of the configuration variables is the motion described by classical mechanics. In the genèric 



34 Isolating constants of the motion correspond to the 'simple' constants of the motion in [60], p.60. 
35 Let us list the introduced dimensions: 

i. - The total dimension is 2(N + 1). 

ii. - The number of the isolating constants is A + 1. 

iii. - The number of the non-isolating constants is N — A. 

iv. - The number of the configuration coordinates is N + 1. 

v. - The dimension of M(po, ■■■,Pa) is 2N — A + 1 
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case, the freqüències of the motions are not rationally dependent (or non-commensurable); thus, the corresponding 
trajectòries fill each torus in a dense way: the motion is ergodic on each torus. 36 Moreover, we can see from eq.(79) 
or in paper [26] that, if the angle-action variables exist, there is a unique equilibrium state on each torus, 

p*{4>) = (2 7r yv+i p ^ ' Pl ' ••■'P Ar ) ( 96 ) 

which is constant in the submanifold M(po, ...,pn); this corresponds to a microcanonical equilibrium on each torus. 



B. Partially non-integrable systems 

In this case, not all the P are isolating constants of the motion; so, there are H and J^. Then, the trajectòries 
must be dense in a domain T>(po, ...,pa) C A4(po, ...,pa) of dimension 2N — A + 1. In fact, if the dimension of 
T>(po, ...,pa) were < 2N — A+1, at least a new global constant would exist and there would be A + 2 global constants; 
but this is impossible since, by its own definition, A + 1 is the total number of these constants. It is quite clear that 
something as a 'thermodynamic equilibrium' classical density must be globally dcfincd in V(po, ...,pa) and, therefore, 
the Ji0 o (0), Jjv-A0 o (0) cannot be explícit variables of p» (</>); 37 then, we must have (as in eq.(79) 38 ) 

= c(gpW> H ^_ HAm P^ H^),...,H Am (97) 
Since T>(po, ...,pa) is contained in the hypersurface defined by 

Hq{4>) = Pa = const., H\{4>) — Pi — const., Ha{4>) = Pa = const. (98) 

C(H ((f>) 7 Hi{4>), H A ((j))) = C(po, P\,...,pa) is the volumc of the portion of the just mentioned hypersurface 
contained in T>(po, ...,pa) (in such a way that p^((f>) would be normalized). Therefore, 

PÏ{4>) = 7^7 " zPÍPo, Pi, -,Pa) = const. (99) 

C(p ,pi, :.,PA) 

and we have found a unique equilibrium for each V(p , ...,pa), namely, for each set of isolating constants of the motion 
or state variables (p n , ...,pa). Now we can repeat the argument of the integrable case. In this case, phase space is 
foliated in submanifolds M{po, ...,pa) of dimension 2N — A + 1. The diffcrcnce is that now not all the coordinates of 
these submanifolds are configuration variables: 39 since there are also momentum variables, in this case we cannot use 
the argument about freqüències not rationally related. Nevertheless, in each T>(po, ...,pa) C À4(po, ■■■tPa) there is a 
unique equilibrium state (99); as a consequence, by using theorem (4.3) of ref. [56] we can conclude that the motion 
is ergodic. 40 

At this point, we may ask ourselves why the density p*(iïo (</>), Hi(<P), Ha{4 } )Ji{4>)i ■••> Jn-a(4>)) looses its J < p i 
variables. The physical reason for this relies on the fact that the space of observables O contains only physical 



36 If the system is unbounded, ergodicity requires that the trajectòries be dense in a domain T>(po, ...,Pn) C M(j>o, ••·,Pjv) (see 
subsection B). 

37 Even if another kind of equilibrium could be defined using local coordinates, it is quite clear that only globally defined 
variables can play the role of thermodynamic variables in eq.(99). 

3S There are not indices 4>i because the variables are the constants of the motion globally defined in a manifold of 2{A + 1) 
dimensions. 

39 Regarding the configuration variables, it is clear that, in the non-integrable case, none of them is a global constant of the 
motion. This fact further reduces the dimension of A4 (or T>). In fact: 

i. - The preserved tori satisfy an irrationality condition ( [58], eq. (3.4.12)); so, the corresponding ratios of the freqüències are 
irrational and the trajectòries are dense in those tori. 

ii. - In the broken tori the trajectòries are chaotic. Nevertheless, if in a particular case there is a configuration variable X that 
turns out to be a global constant of the motion, it can be considered among the 'H'. Then, we will essentially work in the 
manifold X = const. and nothing will change. 

40 Thcorcm (4.3) requires that the evolution be a Frobenius-Perron one. As explained in section V, this is precisely the case. 
In the equilibrium h — ► limit, the motion takes place along classical trajectòries. So, even if p*(<p) is constant in time, the 
partides of the system are in motion. 
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measurable observables. In the integrable case, we have a global ./V + 1 — CSCO whose observables can be globally 
measured in an independent way because they commute. More generally, we can measure only the variables that belong 
to a global set of commuting observables, even if the number of observables is < N+l. In fact, the J could be measured 
in a local N+ 1— CSCO, but they change when they go from one local N+ 1— CSCO to another; therefore, the period 
for making the measurement could be very short (eventually shorter that the period necessary for the measurement 
itself, turning the measurement impossible). In other words, since the classical momenta Ji,^ (</>), Jn-A4>í{<Í>) have 
an ergodic motion, it is reasonable to suppose that the quantum anàlogues Ji^í, Jn-A4>í cannot be really measured, 
even at the quantum level. 41 This means that the set {Ho, Ha} is the relevant measurable global A + 1 — CSCO: it 
is possible to measure only the H, H\, ...,Ha- Then, the isolating constants Ho{4>), ...,Ha{<P) turn out to be the only 
reliable characters in the quantum or the classical play, and the only candidates for thcrmodynamic variables. 

This point can be made in a diffcrent way. In the coordinates , states can only be locally diagonalized. For this 
reason, it is convenient to consider that the unitary operator U of eq.(51) diagonalizes only the indices of the H (that 
we will call r) and does not diagonalize the indices of the J (that we will call m) . Then, we obtain a basis where the 
coordinates of the states read (see (52)) 

Pi^ipirmr'm' = P<pirmm' {oj)S rr i (100) 

With this basis, eqs.(53) and (54) become 

p* = W lim p(t) = V" (L·pé irmm >(u))(ui,rrnrrn'\é i (101) 

irmm' 

and 

Pj= J d ^ P krrn^rm)^{Lürm ! \^ 7=1,2,.., ,4 (102) 

irmm' 

and so on for the rest of the equations. But the m indices can be 'traced away' because the corresponding observables 
cannot be measured, i.c. they cannot be considered classical in a global way. This strategy amounts to consider that 
the J operators and the m indices are inexistent in space O, to the extent that O is the space of all measurable 
observables. This fact is manifested by the following change in the observables of eq.(40) 

O(w)0 ormr ' m ' — ► O(u>),p 0rr 'S mm i (103) 

where we only consider the diagonal term since we are only concerned with the classical part. As we can see, the 
m indices have a 'spherical symmetry' and they measure nothing. Moreover, with variables H and their conjugate 
configuration variables, we can define a 2(A + 1) manifold </>o, and with variables J and their conjugated configuration 
variables, a 2(N — A) manifold <pi rre , in such a way that M — <p Q ® (j)i rre - Then, O(üj),p orr i is globally defined in O 
and 

O s = ^2 j du>0(u)) < f, orr 'd mm '\u)rmr' W)^ = 

irmr'm' 

rr' im rr' 

where 

\ UJrr ')<t>a = X! \ Lürmr ' m )4> i 
im 

is the basis of the new space of observables. Then, the relevant part of eq.(40) reads 



We can only measure dynamical variables when they can be considered as constants in time, at least in the period of 
measurement. If a constant is not global, it is only constant in time in a local coordinate system, i.e. it is not a usual 'physical 
constant'. 
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(P)pl = E J ( ^ UJ P{ iü )<l>irrar'm'0(uj) t j lirmrlml — ^ J dlü pfcrmm' {oj)S rr > 0(üj) <p orr i S n 



E/4e 

r \ im 



P{u)<t>irmm O(u)^ orr (104) 



Calling Eim /'( w )^rmm = p{w)<j, or , namely, making the m— trace of p(w)0 o rrrir'm, we obtain the 'traced' equation 

(Ò)pJ = E / ^(^W 0( w Wr (105) 
where the m indices have disappeared; eq.(lOl) now reads 

P* T = W lim p(i) T = V / (iu>p(u;) 0or (w,r| 0o 

r J 

where T now means 'm— traced' and {(u>,r|0 o } is the dual basis of {|w,r)0 o }. If we work only with the r indices 
(ro, ...,ta), the p of eq.(96) becomes the p of eq.(99) solving the problem. 42 After this coarse-graining that rctains 
only the r variables, the relevant phase space <fio has only 2(A + 1) dimensions. So, we can repeat all the theory with 
only these dimensions and find that, at each manifold M(po,pi, ...,pa), there is only one constant equilibrium state 
PÏ{4>) given by eq.(99) with p(p ,Pi, -,Pa) = p{^)<t> r = Po, r = (pi, ■■■ 1 Pa))- 
The above conclusions are important for two reasons: 

1. - They prové that a system with a continuous evolution spectrum is ergodic in M.(p ai pi, ...,pa) with a constant 
equilibrium density p^ . Therefore, such a system satisfies the two bàsic hypotheses of classical statistical mechanics: 
it is ergodic and microcanonical. 

2. - All the above classical properties can be transferred to the quantum level; therefore, we can define quantum 
ergodic systems as those that satisfy these transferred properties. This definition would be one of the possible 
defimtions of quantum chaos (see [61]): precisely the one which says that a quantum system is chaotic if it has a 
chaotic classical limit ( [45], [47]). This fact is a consequence of our definition of the classical limit. 

C. Non-integrable systems 

Let us now consider the particular case of a non-integrable system where only H((j)) is globally defincd (A = 0). 
Then, we can trace away all the J^: the resulting traced system has a global CSCO {H}. As a consequence, in each 
energy manifold H = lo = const. the motion is ergodic: this is the typical ergodic motion (the drop in the glass of 
wàter, etc). In this case, after tracing away the J ( f >i (i.e. the r-coordinates disappear), the pointer basis only depends 
on H, namely, it is just the eigenbasis of H. 

VII. CONCLUSIONS 

We want to conclude the paper proposing some suggestions for future rescarch. 

i.- We have essentially presented a minimal formalism for quantum chaos, to the extent that our quantum formalism 
satisfies a minimal requirement for such a theory: by definition, a quantum chaotic system has, at least, a classical 
non-integrable system as its classical limit. In fact, this is a necessary but not a sufiïcient condition that any proposcd 
theory of quantum chaos must fulfil. Therefore, our next task is to address the question of whether the set of 
phenomena known under the name of 'quantum chaos'( [45], [47], [51], [55], [61]) can be explained by means of our 
theoretical structure. 



42 At this point, we could say that we are working with a system given by the Ho, H\, Ha and an environment given by the 
J<p i, Jn-A4, - we could consider this situation as the case of an operi quantum system. 
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ii. - Quantum contexts are clearly related with N + 1— CSCOs. We have seen that genèric N + 1— CSCOs are local. 
This might have a relation with well known physical qüestions, as the EPR problem and the Kochen-Speker and Bell 
thcorems (see [52]), where paradoxes arise when we try to describe the quantum system with just one CSCO. 

iii. - In some sense, the equations of quantum physics have a local character ( [62], [63], [64]); we have found that 
this is also the case of the CSCOs: it might be useful to cxplore this analogy. 

iv. - In quantum gravity we have to combine concepts of a local classical theory, general relativity, with thosc of a 
global one, quantum mechanics, and this task has shown to be almost impossible. Now that we have understood the 
local nature of quantum mechanics, the task might become easier. Or rephrased in a more modest way: in quantum 
field theory in curved space time, we have the unsolved vacuum dcfinition problem ( [65], [66], [67]), where the vacuum 
(a clearly local classical general relativity notion) is considered on the global grounds of the quantum field theory. 
Perhaps the present technique might contribute to solve the problem. 
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APPENDIX A: INTEGRABILITY OF THE DIFFERENTIAL EQUATIONS (19) AND (23) 

a.- Let us consider eq.(19), first as a mathematical partial diffcrcntial equation, and then from a physical point of 
view. Calling J^- = A,-, = -B jt eq.(19) reads 

E. dOi ^ dOi „ , . . 

)rf ' 4 ^ + Bi 9ST° (A1) 

where functions Aj, Bj are known. To solve this equation around a hypersurface T> (containing a point fa, that we 
will call (gg°\ ¥>^)) is equivalent to find the trajectòries crossing V such that 



(A2) 



dgp _ _ dq N _ dp q0 _ _ dp qN 
Ao An Bo Bn 

This requires to solve the system of ordinary differential equations 

dqi Ai dq N A N 

dq A dq A Q 

dpqo B dp q N B N 

-r- = ~T = °o,-,-; — = - r = b N A3 
dq a A Q dq A 

where we have taken qo as a parameter. Calling ip = {cp a } = {qi, qN,PqO, ■■·,PqN}, the Lipschitz conditions on 
functions a and b are satisfied if there exists a number M > such that 

2N+1 

\a(q ,<p') -a(q Q ,ip)\ < M ^ \ip' a - ip a \ 

a=l 



2N+1 

\b(q , </) - 6(90, <p)\ < M W a ~ Va\ (A4) 

If functions a and b are continuous and satisfy these conditions in a ncighborhood of a point (q^ , ip^ ) considered 
as the initial conditions, then there exists a unique solution of the system (A3) in a ncighborhood of (p^ and for 
qo £ [<?o°' — c ; q^ + c] for some c > 0. Moreover, if also 
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where the <èj are qo and the 27V + 1 parameters ip a on the surface V (see [46] ) , then we will also find a local solution 
of the partial-differential equation (Al). 

b.- Let us now consider eq.(19) from a physical point of view, in order to see if the Lipschitz conditions are satisfied. 
Introducing the time t — q and, therefore, po = H independent of t, eqs.(A3) read 



dqi 
dt 



dH(q,p) 



dpi 
dt 



dH(q,p) 



dpi ' dt dqi 
namely, the Hamilton cquations. This system is Lipschitzian if there is a M > such that 

dH(Q,P) dH(q,p) 



(A5) 



dp t 



dpi 



< 



8H(Q,P) dH(q,p) 



dqi 



dqi 



M^JÇlQt-qil + lPi- 



(A6) 



in a certain domain V of the phase space. These conditions are fulfilled if the derivatives 9 , d g q í< Qp P ^ ; d df < d p P ^ 

are bounded in the domain V (see e.g. [68] p.141). We can understand the meaning of this requirement in a physical 
example. Let us consider the Hamiltonian 



H(qi 1 q2 1 pi,p 2 ) 



Y + Y + Y + Y + V{quq2 > 



(A7) 



For this particular casc, the derivatives of the potential 9 ^q^ must be bounded in V. Let us consider the potentials 

V(Qi,Q2) = (?i — l2) a - For a > 2, these potentials have their 9 ^q^ bounded in V and the system is Lipschitzian 
(this is the case in quadratic or polynomial potentials used in many models). For a < 2, these potentials have their 
9 dq^dq^ n0 ^ bounded in T> and the system is not Lipschitzian (this is the case of Newton potential (qi — q 2 ) ~ 1 ). But 
realistic physical potentials are always bounded, e.g. the Lenard-Jones potential for finite energies. In fact, there 
is always a central repulsive core. This is not the case of the (mathematical) three bodies problem with Newton 
potentials, provided that all the planets can get infmitely close; but, of course, this is not a physical situation. In 
conclusion, all the physical potentials are Lipschitzian. 

c- Let us now consider the solution of cq.(23) in a powers of h expansion. A simple computation shows that, in 
this case, the problem reduces, for each term, to the one studicd in paragraphs a.- and b.- of this Appcndix, but for 
non-homogencous equations. In fact, 



{Oj(4>), Oj{4>)} m b - {Oi{4>), Oj{4>)} vb + h 2 A^(cb) + ... = o 



Then, if 



we have that 



so, 



OM) = of ] (4>) + tfoy>(4>) + .. 



2n(!)/ 



{oj°V),O< o) (0)} pb + h 2 ({0 ( p(<p),0 { j\ ( f>)} Pb + A (1 H(f>)) + ... = o 



(0), 



{O^ (<(>), 0^m pb + A^(cP) =0,... 

namely, the non-homogeneous version of eq.(Al), which can be solved by mcans of the same method. In this case, 
the domain of the solution will be the intersection of the domains of solution of all these equations. 
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APPENDIX B: SINAI BILLIARD. 



Let us consider the Sinai billiard of fig.l [47]. It is clear that, when the is confined to the inside of the billiard, the 
trajectòries are defined by two independent constants of the motion, H and P x (or H and P y , or P x and P y ), which 
constitute a complete set of local (i.e. in the interior D of the billiard) constants of the motion in involution. When 
the ball strikes the boundaries, it is symmetrically reflected, i.e. the incident angle is equal to the refiected angle, and 
the value of some of the constants of the motion changes: for the two horizontal boundaries, H and P x still constitute 
a completo set of local constants of the motion in involution, but P y changes its sign; for the vertical boundary, H 
and P y still constitute a complete set of local constants of the motion in involution, but P x changes its sign. 

Without modifying the physical characterization of the example, we can replace the rigid walls with infinitcly 
hcight potential barriers of width d, namely, the potentials V(x),V(y) and V(r) of figuro 2. (e.g. V(x) bchaves as 
V(0) = 0, V(0) = 0, V(—d) — ► oo). Due to the symmetry of the potentials (translation symmetry for V(x) and V(y), 
rotation symmetry for V(rj), the reflections are still symmetric, i.e. the ball climbs the potential walls and then falls 
down with symmetrical motion. Calling D\ and D 3 the domains in the potential of the x walls, D 2 that of the y wall, 
and £>4 that of the curved wall, we see that x is a eyelic variable in D\ and D3, y is a eyelic variable in D2, and 6 is 
a eyelic angular variable in D4 . Thcrcfore, we have the following local constants of the motion in cach domain: 



Do 


H 


P x 


D 1 


H 


P x 


D 2 


H 


Py 


D 3 


H 


P x 


Di 


H 


Pe 



In summary, we have found five domains, each one with two constants of the motion in involution. If d — ► 0, we go 
from fig.2 to fig.l. 
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Figure Caption: 

fig.l A Sinai billiard. 

fig.2 A Sinai billiard with potential barriers. 
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